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Some of the pitfalls of automatic computation are described as well as their possible avoidance. 
A detailed description is given of the methods for computing the error and complementary error 
function. The implementation of the method is given in a double precision ASA FORTRAN program. 
The test ("driver") program and test results are also given. 
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1. Introduction 

AMS55[1] 1 — the Handbook of Mathematical Functions — in addition to defining functions 
in terms of their integrals, power series, asymptotic expansions, etc., also supplies tables and 
approximating polynomials as well as methods of computing values both within and outside the 
tabular range. With the present trend to large scale computer systems, questions frequently 
arise concerning methods of computing the tabular values themselves. In many, if not most cases, 
the method used for computing the table is certainly not the best method for computing discrete 
values. While many specialized methods and programs exist [3, 6, 7, 8|, as a supplement to the Hand- 
book and to satisfy general needs, a systematic effort is being made to supply methods of computing 
special functions, accurate within machine precision, valid throughout the entire range of the func- 
tion, efficient, and readily programmable and adaptable for varying precision and computers. Parti 
of the present paper is devoted to the error, probability, and related functions. 

2. General Methods 

Methods used in computing functions vary from expansion in series — power, Taylor, asymp- 
totic, etc.; integration of the differential equation; numerical quadrature; evaluation of continued 
fractions; polynomial or rational approximations; recurrence techniques; iteration methods; stored 
tables with the use of interpolation; addition or duplication formulas, etc. Numerical quadrature 
and integration of the differential equation, while ideal for preparing a table, are generally in- 
efficient for obtaining discrete values. Since the coefficients in approximating polynomials vary 
with precision, this method will be avoided if possible. The methods chosen for computing the func- 
tions will be those which avoid some of the common pitfalls or those pitfalls not easily remedied. 

3. Common Pitfalls and Their Avoidance 

Computation with alternating series leads to a loss of significant figures due to either or 
both subtraction and accumulated round-off errors. The use of Taylor series leads to difficulties 
of premature truncation if any of the derivatives possess zeros. The use of slowly convergent 
series leads to the difficulty of a heavy accumulation of round-off errors as well as of devising a 
satisfactory cut-off procedure. In computing with such series one might impose a condition on 



1 Figures in brackets indicate the literature references on p. 212. 
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the number of terms to be used in the computation. However, this number varies with precision 
and the method would impose a timing penalty in certain ranges. After the nth term has satis- 
fied the convergence criteria, the sum of the nth to the (n-\-p)th term should be checked for its 
contribution. The loss of significant figures away from or at the zeros of the functions will be 
noted by comparing the largest term with the sum and the convergent term. This method of cur- 
tailing the series will lead to time savings and frequently avoid overflow and underflow conditions. 
In using asymptotic expansions the test for convergence should precede the test for divergence. 

In computing with continued fractions [2] to avoid premature truncation and or possible 
overflow or underflow difficulties, the convergents should only be obtained after the numerator 
and denominator of the ratio are tested for zeros. In this circumstance, the test for convergence 
should likewise be omitted. In continued fractions where the partial numerators and denominators 
are always positive, no zeros arise and the value of the continued fraction always lies between 
two consecutive convergents, except possibly for round-off errors in the calculation. Use of the 
transformed continued fraction containing only the even or odd convergents of the original frac- 
tion, while more efficient generally, may lead in certain cases to a loss of significant figures and to 
overflow conditions more quickly. 

In recurrence techniques, there may be a loss of significant figures or a heavy buildup of 
round-off errors. Forward and backward techniques or a combination of both may be useful [4, 
5]. Attention must be paid to the proper choice of normalization factor, avoiding the use of values 
at or near zeros. 

In using iterative methods, common pitfalls are the imposition of too stringent a convergence 
criteria and/or a failure to impose a limit on the maximum number of iterations. 

4. General Procedures 

In addition to detailing the methods of computing, the ASA FORTRAN programs will be 
presented, prepared in such a way that only a minimal number of changes need be made to ac- 
commodate the change in precision. Initially the test ("driver") program will also be given as well 
as test results and timing in the hope of detecting software and hardware errors. 

The authors, hoping to satisfy the needs of all users, earnestly solicit suggestions for improve- 
ments and comments relative to testing on other computers. In this way it is hoped that an ex- 
cellent library of special functions will be established. 
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Automatic Computing Methods for Special Functions 

Part I. Error, Probability, and Related Functions 

1. Definitions 

Error Function 

erf x = -^ ( X e' t2 dt erfc * = ~7= | e~ t2 dt 

V77 J {) \Ztt J* 

= 1 — erf x 

erf (— x) = — erf x erfc (— x) = 2 — erfc* 

Normal or Gaussian Probability Function 

V2tt V27T J-* •>-•*• 

= erf (x/V2) 

V27T J-* V27T J* 

= P Z(t)«fc =\~Z(t)dt 

=i[l + erf U/V2)] = 1-P(%)=ierfc(*/V2) 

Confluent Hypergeometric Function 
M(if,-x 2 ) = (Vif/2^) erf* 

M(l, f , * 2 ) = (V^I2x)e** erf a 

£/(i,4,* 2 ) = V^erfc* 
Incomplete Gamma Function 

r(i, a 2 ) =2 J ~e- (2 dt= V^erfc* 

2. Method 

x, a real variable 

NBC = the number of binary digits in the characteristic of a floating point number 
NBM = accuracy desired or the maximum number of binary digits in the mantissa of a floating 
point number 
For x = 0, erf x = 0; erfc x = 1.0. 

213 



The power series is used for \x\ =£ ULPS(= 1), the upper limit for the power series 

2x 



V^ e ' &1-3 . . . (2/V+l) 



77 



2^+i 



?r 2 PWR 

where r =l, D = 1, ^v=~y T s ^ = ^^ T N - i; D N =D N - t + 2; TV ^ 1. 



Since min V T N = 1, it is valid to obtain M for the final term by testing for 



ru^T0LER(=2 NBM ). 

NOTE: In the program, in order to minimize round-off errors, summing starts at N=l. 

One then obtains erfc x— 1 — erf x, erf (— x)=— erf x, erfc (— x)=2— erfc x. 
The continued fraction 



erfc 



e-* 2 T 2* 



1-2 



3-4 



V^|_2* 2 +l-2* 2 + 5-2x 2 + 9- 



-X* 



Vtt' 



is used for 1 < |x| ^ ULCF (= 0. 83 [2 (N BC ~ x )/2 ]), the upper limit (machine bound) for the continued 
fraction. The above continued fraction yields only even convergents, forming a monotonically 
increasing sequence, of the following continued fraction 



erfc x- 



1 1/2 1 3/2 

x + x~\- x + x+ . 



Setting 

F-i = 0, F = 2*, G_, = 1, G =2* 2 + 1, r, = 1, B, = G + 4 with r. v = r V -i + 2, B N =B N -i + 4, iV>l, 
one obtains 



A n =-Wn(Wn+1),F n =AhF n -* + BhFn- u Gh=A&m + Bi& n -uN*1. 

F N /G N 



If 



1 



Fn-iIGn-i 



TOLER(=2-nbm) 9 then Fn j Gn = F. If F N -ilG N -i denoted by PREV^F. V /G* 



the sequence is permitted to continue after first checking for boundedness (G A < scaling factor, 
SCF = 2 C2 - NBC where C2 = 2 NBCl ). If PREV > F N IG N , PREV = F. One then obtains the quantities 

erf x= 1 — erfc x, erf ( — x) =— erf x, erfc ( — %) =2 — erfc %. 

For |jc | > 0.83 [2 (NBC_1) / 2 ], erf |x| = 1, erfc \x\ =0, the limiting values of the functions. 

Beyond \x\ = 1, a loss of significant figures will occur in subtraction in obtaining erfc x, and 
in addition round-off errors begin to accumulate more rapidly. If only erf x is needed, a time saving 
is achieved by extending the power series to \x\ =2. The continued fraction should then be used 
until the appropriate critical value is reached. 

While the continued fraction is convergent for x > 0, the number of convergents needed in- 
creases rapidly and without scaling underflow or overflow occur for x nearing zero. For NBC = 8, 
NBM = 27, underflow occurred near x = 0.91. The partial numerators and denominators F\-i, 
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Gjv_i, F\- 2 and Gs->i are scaled but greater time saving is achieved by using the power series. If 
the continued fraction for \x\ < 0.61 both F\ and G\ must be tested for boundedness. 
maximum precision attainable with the use of the asymptotic expansion 



one uses 
The 



srfc x = 



1 , ^ (-l)- v l-3. . . (27V- 1) 



(2x 2 )' 



for 



The 
of NBM 



# = 4.6, |e| ^ 1 unit in the 9th significant figure, 

# = 6.7, |e| ^ 1 unit in the 19th significant figure. 

following table gives an indication of the number of terms needed for particular values 





Maximum number of terms 




NBM = 27 


NBM = 60 


Power series — erf 1 


13 

24 
9 
6 
3 
3 


21 


Continued fraction — erfc 1 


124 


Asymptotic expansion — erfc 6.7 


47 


erfc 9.3 

Continued fraction — erfc 6.7 


17 
8 


erfc 9.3 


6 







3. Range 

The range covered is limited by the machine bound for the exponential function. 
For |*| ss 0.83 [2 (NH(,)/2 ], the computations are performed with \x\^ 9.3 for NBC = 8, |*|^ 26.5 
for NBC = 11. 

For |*| > .83[2 (NBC - 1)/2 ],erf |*| = 1, erfc |*|=0. 

4. Accuracy 

At least (NBM-I-3) significant binary digits where I is the number of binary digits repre- 
senting the integer part of x 2 . (This is essentially the accuracy of the exponential routine, with 
round-off error allowance.) 

5. Precision 

Variable — by setting the desired NBM. 

6. TIMING (Seconds- Univac 1108) 



For the range 




NBM = 27 


NBM = 60 


Method 


^ = 0(.005)1 


(201 values) 


0.08 
.16 


0.22 
.531 
.13/ 


Power series 


*=1 (.05)9.4 


(169 values) 


Continued 


* = 9.2(.2)26.6 


(88 values) 


fraction 








Maximum Time/Eva 




.003 


.015 
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388-513 O - 70 - 7 



7. Testing 

To avoid errors due to incorrect arguments, the beginning argument, the interval, and end 
argument were read in as scaled integers. The generation was performed with integer values, 
double precision conversion made, then single precision values obtained. 

The double precision results were checked, when possible, against those obtained by using 
the polynomial or rational approximations of Clenshaw [l], 1 Cody [2], Hart [3], and Strecok [6]. In 
addition, special values were checked by the asymptotic expansion and numerical integration. 
The single precision results were then checked against the double precision results. In all cases 
the results obtained agreed within the reported accuracy. 

For those desiring both functions accurately with limited precision throughout the range, 
the approximations of Cody were found the most efficient and accurate. 

8. Critical Values 

erf c 9.3-0.16532 44184 030157 (-38) (NBC = 8) 
erfc 26.5-0.22109 07664 26373 (-306) (NBC=11) 



erf 4.33 
erf 4.82 
erf 5.68 
erf 6.43 



= 1 



|e|<l(-9) NBM-27 
\e\ <1(-11) NBM = 36 
|e|<l(-15) NBM = 48 
|c| <1(-19) NBM-60 



9. Special Constants 

2(tt)"£ = 1.12837 91670 95512 57389 61589 03121 

(tt)*/2 = 0.88622 69254 52758 01364 90837 41670 

(fr)"* =0.56418 95835 47756 28694 80794 51560 

tt* = 1.77245 38509 05516 02729 81674 83341 

(2tt)"* =0.39894 22804 01432 67793 99460 59934 

(2tt)* =2.50662 82746 31000 50241 57652 84811 

77=3.14159 26535 89793 23846 26433 83279 

(2)~* =0.70710 67811 86547 52440 08443 62104 

2* = 1.41421 35623 73095 04880 16887 24209 

erf 1=0.84270 07929 49714 86934 12206 
erf p = 0.5 

p= 0.47693 62762 04469 87338 



1 Figures in brackets indicate the literature references on p. 217. 
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2" 27 - 0.74505 80596 92382 8125 ( - 8) 

2 "=0.14551 91522 83668 51806 64062 5 (-10) 

2" 48 -0.35527 13678 80050 09293 55621 33789 0625 (-14) 

2™ =0.86736 17379 88403 54720 59622 40695 95336 91406 25 (-18) 

10. References 
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5, London, 1962. 
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f 3 J Hart. John F. et al., Computer Approximations, The SIAM Series in Applied Mathematics (John Wiley & Sons, Inc., 

New York, 1968). 
[4| National Bureau of Standards, Tables of Normal Probability Functions, Applied Mathematics Series 23 (U.S. Government 

Printing Office, Washington, D.C., 1953). 
[51 National Bureau of Standards, Tables of the Error Function and its Derivative, Applied Mathematics Series 41 (U.S. 

Government Printing Office, Washington, D.C., 1954). 
[6| Strecol , Anthony J., On the calculation of the inverse of the error function, Math, of Comp. 22, 144-158 (1968). 
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1* SUBROUTINE ERRINT (X, ERF, ERFC ) 

2* C (ASA FORTRAN) 

3* C USAGE. CALL ERRINT (X, ERF, ERFC) 

4* C ENTERING WITH X(A REAL VARIABLE) , ONE OBTAINS BOTH ERF X 

5* C AND ERFC X TO THE MAXIMUM MACHINE ACCURACY POSSIBLE. 

6* C 

7* C THE PROGRAM, AS SET UP BELOW 

8* C ( 1 ) WITH THE DOUBLE PRECISION TYPE STATEMENT , 

9* C ( 2 ) NBC=11 , NBM=60 ON THE FIRST DATA CARD AND 

10* C (3) THE REMAINING CONSTANTS ON THE SECOND AND THIRD DATA 

11* C CARDS GIVEN AS DOUBLE PRECISION CONSTANTS, YIELDS DOUBLE 

12* C PRECISION RESULTS ON THE UNIVAC 1108 SYSTEM. 

13* C 

14* C SINGLE PRECISION RESULTS ARE OBTAINED BY 

15* C ( 1 ) DELETING THE TYPE STATEMENT 

16* C ( 2 ) SETTING NBC=8 , NBM=27 ON THE FIRST DATA CARD 

17* C ( 3 ) CHANGING THE D ' S TO E ' S ON THE SECOND AND THIRD DATA 

18* C CARDS. 

19* C (4) CHANGING FUNCTION TYPE-DABS TO ABS , DEXP TO EXP 

20* C FOR OTHER VALUES OF NBM, THE CONSTANT ON THE THIRD DATA 

21* C CARD SHOULD BE GIVEN TO THE CORRESPONDING PRECISION . 

22* C 

23* C CAUTION. THE SUBROUTINE CANNOT READILY BE ADAPTED TO 

24* C COMPUTE THE ERROR FUNCTION FOR A COMPLEX 

25* C ARGUMENT. 

26* C ERF(X) =(2/SQRT(PI) ) *INTEGRAL FROM TO X OF EXP (— T**2)DT 

27* C ERFC ( X ) = ( 2/SQRT ( PI ) ) *INTEGRAL FROM X TO INFINITY OF 

28* C sxP (— T**2) DT=1— ERF (X) 

29* C 

30* C ERF(-X) =-ERF(X) 

31* C ERFC(— X)=2— ERFC(X) 

32* C 

33* C NBC=NO. OF BINARY DIGITS IN THE CHARACTERISTIC OF A FLOATING POINT 

34* C NUMBER. 

35* C NBM=ACCURACY DESIRED OR MAX . NO . OF BINARY DIGITS IN THE MANTISSA OF 

36* C A FLOATING POINT NUMBER . 

37* C 

38* C METHOD . POWER SERIES FOR ABS ( X ) . LE . 1 , ULPS , UPPER LIMIT FOR 

39* C SERIES , CONTINUED FRACTION FOR ABS ( X ) . GT . 1 AND . LE . 

40* C ULCF (UPPER LIMIT OF CONTINUED FRACTION) 

41* C 

42* C RANGE. ABS (X) . LE . ULCF ULCF= . 83* ( 2 . ** ( (NBC-1 )/2 ) ) 

43* C ABS ( ULCF ) APPROX . 9.3, NBC=8 

44* C 26.5,NBC=11 

45* C BEYOND THIS RANGE THE LIMITING VALUES OF THE FUNCTIONS 

46* C ARE SUPPLIED ERF(INF)=1 . 0, ERFC ( INF )=0. 

47* C 

48* C ACCURACY . ROUTINE WILL RETURN ( NBM-I-3 ) SIGNIFICANT BINARY DIGITS 

49* C WHERE I IS THE NUMBER OF BINARY DIGITS REPRESENTING THE 
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50* 
51* 
52* 
53* 
54* 
55* 
56* 
57* 
58* 
59* 
60* 
61* 
62* 
63* 
64* 
65* 
66* 
67* 
68* 
69* 
70* 
71* 
72* 
73* 
74* 
75* 
76* 
77* 
78* 
79* 
80* 
81* 
82* 
83* 
84* 
85* 
86* 
87* 
88* 
89* 
90* 
91* 
92* 
93* 
94* 
95* 
96* 
97* 
98* 
99* 
100* 
101* 



C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 



c 
c 
c 



c 
c 



c 
c 
c 



c 
c 
c 

c 
c 
c 



PRECISION. 

MAXIMUM TIME. 
(SECONDS) 



INTEGER PART OF X**2. (THIS IS ESSENTIALLY THE ACCURACY 
OF THE EXPONENTIAL ROUTINE. ) 

VARIABLE-BY SETTING DESIRED NBM . 

UNIVAC 1108-EXEC II 

NBM=27 NBM=60 
.003 .015 



11 



DOUBLE PRECISION AN, BN, CONS, CI, DN, ERF, ERFC, F, FN, FNM1 , 

1 FNM2, FOUR, GN,GNM1 GNM2 , ONE, PREV, PWR, RNBC, SCF, SUM, 

2 TN, TOLER, TRRTPI , TWO, ULCF, ULPS , WN , X, Y, YSQ 

DATA NBC, NBM/11, 60/ 

DATA ONE, TWO, FOUR, ULPS, CONS/1. DO , 2 . DO, 4 . DO, 1 . DO, . 83D0/ 

DATA TRRTPI/1 . 128379167095512574D0/ 

RNBC=NBC 
T0LER=TW0** (-NBM) 

TEST ON ZERO 

IF (X) 2,1,2 

ERF=0 
ERFC=ONE 

RETURN 

Y=DABS(X) 
YSQ=Y ** 2 

IF (Y-ULPS) 3,3,4 

MAXIMUM ARGUMENT 

C1=TW0 ** ( (RNBC-ONE)/TWO) 
ULCF=CONS * CI 

SCALE FACTOR 
SCF=TWO** (CI** 2— RNBC) 

LIMITING VALUE 

IF (Y-ULCF) 10,10,11 

ERF=ONE 
ERFC=0 
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102* 




103* 


C 


104* 


c 


105* 


c 


106* 




107* 




108* 




109* 




110* 




111* 




112* 




113* 


c 


114* 


c 


115* 


c 


116* 




117* 


c 


118* 




119* 




120* 


c 


121* 


c 


122* 


c 


123* 




124* 


c 


125* 




126* 




127* 


c 


128* 




129* 


c 


130* 


c 


131* 


c 


132* 




133* 




134* 




135* 




136* 


c 


137* 




138* 


c 


139* 




140* 




141* 


c 


142* 




143* 


c 


144* 




145* 




146* 


c 


147* 




148* 


c 


149* 


c 


150* 


c 


151* 




152* 




153* 


c 



GO TO 7 

METHOD -POWER SERIES 

3 SUM=0 
DN=ONE 
TN=ONE 
PWR=TWO*YSQ 

6 DN=DN + TWO 
TN=PWR*TN/DN 
SUM=TN + SUM 

TOLERANCE CHECK 

IF (TN-TOLER) 5,6,6 

5 ERF-^ ( SUM 4- ONE ) *TRRTPI *Y*DEXP (— YSQ ] 
ERFC=ONE — ERF 

NEGATIVE ARGUMENT 

7 IF (X) 8, 9, 9 



ERF-- 


ERF 




ERFC- 


TWO- 


-ERFC 


RETUR1N 








METHOD ( 



CONTINUED FRACTION 

10 FNM2 = 

GNM2 = 0NE 

FNM1-TW0 * Y 

GNM1 = TW0 * YSQ + ONE' 

PREV = FNM1/GNM1 

WN-ONE 

BN = GNM1 4- FOUR 

14 AN--WN* (WN+ONE) 

FN-BN * FNM1 + AN * FNM2 
GN-BN * GNMl-f AN * GNM2 

F-FN/GN 

TOLERANCE CHECK 

IF (DABS (ONE- (F/PREV) )— TOLER) 12, 12, 13 
13 IF (PREV-F) 17, 17, 18 

BOTH FN AND GN MUST BE TESTED IF ABS (X) .LT. .61 
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154* 




17 


IF (GN .LT. SCI 


155* 


C 






156* 


C 




S 


157* 


C 






158* 




15 


FN-FN/SCF 


159* 






GN=GN/SCF 


160* 






FNM1-FNM1/SCF 


161* 






GNM1=GNM1/SCF 


162* 


C 






163* 




16 


FNM2 = FNM1 


164* 






GNM2-GNM1 


165* 






FNM1 = FN 


166* 






GNM1 = GN 


167* 


C 






168* 






WN-WN-fTWO 


169* 






BN=BN + FOUR 


170* 






PREV-F 


171* 






GO TO 14 


172* 


C 






173* 




18 


F-PREV 


174* 




12 


ERFC-F*DEXP (- 


175* 






ERF = ONE-ERFC 


176* 


C 






177* 






GO TO 7 


.178* 


C 






179* 






END 



SCALING 



-YSQ)*TRRTPI/TWO 
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1* 

2* 
3* 
4* 
5* 
6* 
7* 
8* 
9* 
10* 
11* 
12* 
13* 
14* 
15* 
16* 
17* 
18* 
19* 
20* 
21* 
22* 
23* 
24* 
25* 
26* 
27* 
28* 
29* 
30* 
31* 
32* 
33* 
34* 
35* 
36* 
37* 



C 
C 
C 
C 
C 
C 
C 



C 
C 



TEST PROGRAM FOR DOUBLE PRECISION ERROR FUNCTION AND COMPLI- 
MENTARY ERROR FUNCTION (UNIVAC 1108) 

CALL ERRINT (X, ERF, ERFC) 



DOUBLE PRECISION DN , DSCALE , ERF, ERFC, X 

DIMENSION M(12) 

DIMENSION X (500), ERF (500), ERFC (500) 

DATAM(l), M(2), M(3), M(4), M(5), M(6), M(7), M(8), M(9J 

1 M(10) , M(ll) , M(12)/0, 5, 100, 100, 10, 5, 95, 10, 

2 95, 5, 270, 10/ 
J = l 

N = M(J) 

NINCR=M(J + 1) 
NENDX=M(J + 2) 
DSCALE = M(J + 3) 

1 = 1 

DN=N 

X(I) =DN/DSCALE 

IF (N-NENDX) 2, 3, 3 

1 = 1 + 1 

N = N + NINCR 

GO TO 5 



3 DO 4 L=l, I 

4 CALL ERRINT (X 

WRITE (6, 1002) 
DO 6 L=l, I 



ERF (L) , ERFC (L] 



38* 




6 


WRITE (6, 1001) X (L), ERF (L), ERFC (L) 




39* 


C 








40* 






IF (J .EQ. 9) STOP 




41* 






J = J + 4 




42* 






GO TO 1 




43* 


C 








44* 




1001 


FORMAT (3D26. 18) 




45* 




1002 


FORMAT (1H1, 11X, 1HX, 23X, 6HERF(X), 20X, 


7HERFC(X) ,/) 


46* 


C 








47* 






END 
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ERF(X) 



ERFC(X) 



500000000000000000-001 
1O0000000000000000 I 000 
150000000000000000 + 000 
200000000000000000 + 000 
250000000000000000 + 000 
+ 000 
+ 000 



450000000000000000 + 000 
500000000000000000 + 000 
549999999999999999 + 000 
599999999999999999 + 000 
+ 000 
+ 000 

750000000000000000 + 000 
799999999999999999 + 000 
849999999999999999 + 000 



+ 001 



.563719777970166236-001 
.112462916018284892 + 000 
.167995971427363495 + 000 
.222702589210478454 + 000 
.276326390168236935 + 000 
.328626759459127429 + 000 
.379382053562310322 + 000 
.428392355046668455 + 000 
.475481719786923677 + 000 
.520499877813046541 + 000 
.563323366325108957 + 000 
.603856090847925924 + 000 
.642029327355671843 + 000 
.677801193837418476 + 000 
.711155633653515133 + 000 
.742100964707660485 + 000 
.770668057608352530 + 000 
.796908212422832131 + 000 
.820890807273277944 + 000 
.842700792949714871 + 000 



i+OOl 

.943628022202983377 \ 000 
.887537083981715109 + 000 
.832004028572636506 + 000 
.777297410789521546 + 000 
.723673609831763065 + 000 
.671373240540872572 + 000 
.620617946437689678 + 000 
.571607644953331545 + 000 
.524518280213076324 + 000 
.479500122186953459 + 000 
.436676633674891044 + 000 
.396143909152074076 + 000 
.357970672644328157 + 000 
.322198806162581524 + 000 
.288844366346484867 + 000 
.257899035292339516 + 000 
.229331942391647470 + 000 
.203091787577167870 + 000 
.179109192726722056 + 000 
.157299207050285130 + 000 



ERF(X) 



ERFC(X) 



+ 001 

150000000000000000 + 001 
+ 001 
+ 001 

300000000000000000 + 001 
350000000000000000 + 001 
+ 001 
+ 001 
+ 001 

549999999999999999 + 001 
599999999999999999 + 001 
+ 001 
+ 001 
+ 001 

799999999999999999 + 001 
849999999999999999 + 001 
+ 001 
+ 001 



.842700792949714871 + 000 
.966105146475310728 + 000 
.995322265018952735 + 000 
.999593047982555043 + 000 
.999977909503001415 + 000 
.999999256901627659 + 000 
.999999984582742101 + 000 
.999999999803383957 + 000 
.999999999998462541 + 000 
.999999999999992643 + 000 
.999999999999999979 + 000 
1 + 001 
1 + 001 

.100000000000000000 + 001 
1 + 001 
1 + 001 

.100000000000000000 + 001 
.100000000000000000 + 001 



.157299207050285130+000 
.338948535246892735-001 
.467773498104726582-002 
.406952017444958938-003 
.220904969985854420-004 
.743098372341412774-006 
.154172579002800195-007 
.196616044154288756-009 
.153745979442803493-011 
.735784791797439867-014 
.215197367124989151-016 
.384214832712064783-019 
.418382560777941505-022 
.277664938603056939-025 
.112242971729829286-028 
.276232407133377193-032 
.413703174651381081-036 
.376921448565488099-040 
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ERF(X) 



ERFC(X) 



+ 001 
+ 002 

105000000000000000 + 002 
110000000000000000 + 002 
115000000000000000 + 002 
120000000000000000 + 002 
125000000000000000 + 002 
130000000000000000 + 002 
135000000000000000 + 002 
140000000000000000 + 002 
145000000000000000 + 002 
150000000000000000 + 002 
155000000000000000 + 002 
160000000000000000 + 002 
165000000000000000 + 002 
170000000000000000 + 002 
175000000000000000 + 002 
180000000000000000 + 002 
185000000000000000 + 002 
190000000000000000 + 002 
195000000000000000 + 002 
200000000000000000 + 002 
205000000000000000 + 002 
210000000000000000 + 002 
215000000000000000 + 002 
220000000000000000 + 002 
225000000000000000 + 002 
230000000000000000+002 
235000000000000000 + 002 
240000000000000000 + 002 
245000000000000000 + 002 
250000000000000000 + 002 
255000000000000000 + 002 
260000000000000000 + 002 
264999999999999999 + 002 
270000000000000000 + 002 



+ 001 
+ 001 

100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
100000000000000000 + 001 
+ 001 
+ 001 
100000000000000000 + 001 

4 001 

+ 001 

100000000000000000 + 001 
1000000000000000004 001 
100000000000000000 4 001 
+ 001 
+ 001 

100000000000000000 + 001 
100000000000000000 + 001 
1000000000000000004 001 
1000000000000000004001 
1000000000000000004001 
100000000000000000+001 
100000000000000000 + 001 
100000000000000000 + 001 
+ 001 
+ 001 



.376921448565488099-040 
.208848758376254514-044 
.703592809017752444-049 
.144086613794369515-053 
.179330964357678251-058 
.135626116920590467-063 
.623194278197991275-069 
.173955731546672518-074 
.294943311325798927-080 
.303722984775031350-086 
.189939594197950386-- 092 
.721299417245120960-099 
.166320164004887325-105 
.232848575157153193-112 
.197913057555326932-119 
.102122801509426159-126 
.319886381234348273-134 
.608236923181640276-142 
.701996157498568456-150 
.491772283925648080-158 
.209095414792273150-166 
.539586561160790480-175 
.845084236957246622-184 
.803245387102246359-193 
.463333653965844984-202 
.162190586093347461-211 
.344534886046460488-221 
.444126594808806166 - 231 
.347405949564997476-241 
.164898258315193574-251 
.474936126406738376-262 
.830017257119653210-273 
.880166269072796021 - 284 
.566319240885615157-295 
.221090766426373696-306 
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